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We demonstrate chaotic behavior of timelike, nuU and spacehke geodesies in non- homogeneous vacuum pp-wave solutions. 
This seems to be the first known example of a chaotic motion in exact radiative spacetime. 
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In the context of general relativity the first system 
for which a chaotic behavior of solutions to the Einstein 
equations has been recognized and thoroughly studied 
were Bianchi IX cosmological models (see for example 
101, 1^, and references therein). Complicated nonlinear 
effects also occur in systems with coupled gravitational 
and scalar fields |^ - (g. 

Another type of problems providing nonlinear dy- 
namical systems in general relativity are the studies 
of geodesic motion in given spacetimes. In particular, 
chaotic behavior of geodesies in the relativistic analogue 
of the two fixed-centres problem (modelled by extreme 
black holes) was examined in ||^ - . Chaotic geodesic 
motion was also found in Schwarzschild spacetime 
- [0, in some static axisymmetric spacetimes \lt] 
and in topologically non-trivial Robertson- Walker uni- 
verse 0, ^ 

Here we investigate motion in exact gravitational 
waves, namely in the widely known class of vacuum pp- 
waves [ p^ , the metric of which can be written in the 
form 



ds^ = 2 dCdC - 2 dudv -if + f) di 



(1) 



where f{u, ^) is an arbitrary function of u and the 
complex coordinate ^ spanning the plane wave surfaces 
u — const. When / is linear in C, the metric (|l|) 
represents Minkowski universe. The case / = g{u)C'^ 
describes plane gravitational waves ('homogeneous' pp- 
waves) which have thoroughly been investigated (see 
for Refs). This simple example of an exact radiative 
spacetime has also been used for construction of sand- 
wich and impulsive waves Q. 

However, here we wish to study geodesies in more 
general, non-homogeneous vacuum pp-waves and demon- 
strate their chaotic behavior. The geodesic equations for 
(HI) are 



ii ^ U ~ const. , 



, 



(2) 
(3) 
(4) 



where dot denotes d/dr with r being an affine parameter. 
Assuming also a condition normalizing the four-velocity 



such that UfjU'^^ = e where e ~ — 1,0,+1 for timelike, 
null or spacelike geodesies, we get 



2CC -(,/ + ./)[/'-£ 



(5) 



We consider [/ ^ (for f/ = Eqs. (|)-(|) can be 
integrated yielding only some trivial geodesies). By dif- 
ferentiating Eq. (^ and using (H) we immediately obtain 
(^ which can thus be omitted. Hence it suffices to find 
solutions of Eq. (H) since v{t) can then be obtained by 
integrating Eq. (||), and u{t) — Ut + uq. 

The remaining equation ^ has the same form for 
timelike, null and spacelike geodesies. Introducing real 
coordinates x and y by C = a; -I- iy we get a system which, 
for / independent of u, follows from the Hamiltonian 



H=^{pl+pl)+V{x,y) 



(6) 



where the potential is V{x,y) — ^U^TZef. For non- 
homogeneous pp-wave spacetimes given by / = C^", C = 
const. > 0, n = 3, 4, • • •, the corresponding potential 



V{x,y) = ^CU^TZeC 



(7) 



is called 'n-saddle'. It can by visualized in polar coor- 
dinates p, (j) where C = pexp(i0), in which it takes the 
form V{p, cj)) = \CU^ cos(n0). 

Now, it was shown previously by Rod, Churchill and 
Pecelli in a series of papers j2^] - that motion in the 
Hamiltonian system (|^) with the n-saddle potential (|^) 
is chaotic. 

Let us first briefly summarize their results for the sim- 
plest case n — i. The corresponding potential (after 
removing an unimportant multiplicative factor by a suit- 
able rescaling of r) 



V{x,y)^\x^ 



xy 



(8) 



is called a 'monkey saddle'. Interestingly, this is a spe- 
cial case of famous Henon-Heiles Hamiltonian which 
is known to be a 'textbook' example of a chaotic sys- 
tem (however, their quadratic terms are absent in our 
case). This particular case of the Henon-Heiles Hamil- 
tonian has been investigated by Rod ||2^. He con- 
centrated on bounded orbits in the energy manifolds 
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H{x,y,Px,Py) = E > 0. The homogeneity of V guaran- 
tees that the orbit structure for any two positive values 
of E is isomorphic modulo a constant scale factor and 
adjustment of time: x^x — Xx,y-^y — Xy and 
T ^ T — t/\/A results in E E = X'^E. Therefore, 
without loss of generality one can restrict to one partic- 
ular value of E. 

In order to describe the topological structure of all 
bounded orbits Rod first constructed three basic unstable 
periodic orbits (denoted by Ilj) which are isolated invari- 
ant sets for the flow. The region in which these bounded 
orbits occur can be decomposed into three disjoint cells 
Rj (see or for details); each contain only one 
orbit Ilj and no other bounded orbits. Hence, Ilj is the 
only invariant set in Rj, i.e. it is isolated. 

Subsequently, Rod investigated orbits asymptotic to 
basic orbits Ilj as r ^ ±c» and showed that these asymp- 
totic sets intersect transversely. This gives the existence 
of orbits that 'connect' the orbits Ilj: they are homo- 
clinic (asymptotic to the same periodic orbit in both 
time directions) or heteroclinic (asymptotic to two dif- 
ferent periodic orbits, one in each time direction). It is 
the existence of these orbits that illustrates complicated 
chaotic structure of the flow. 

The topology of possible orbits in phase space can also 
be represented by symbolic dynamics given here by a set 
of biinfinite sequences, s = ■ ■ ■ , Sk, Sfc+i, Sk+2, ■ ■ ■, where 
Sk S {1,2,3}, Sk Sk+i- Using a topological version of 
the Smale horseshoe map, it was shown in p2[ that to any 
biinfinite sequence of symbols {1, 2, 3} there exists an un- 
countable number of bounded orbits running through the 
blocks Rj in the prescribed order as r goes from — oo to 
-\-oo . Also, the flow admits at least a countable number 
of non-degenerate homoclinic and heteroclinic orbits. 

Rod remarked that these results could be refined if the 
unstable periodic orbits Ilj were known to be hyperbolic 
so that they would admit stable and unstable asymp- 
totic manifolds. Consequently, to each periodic symbol 
sequence there would correspond a countable collection 
of periodic orbits. The difficulties in proving the hypcr- 
bolicity of Ilj were subsequently overcome in . In \ , 
summarizing and generalizing some previous results P5 |, 
the Hamiltonian (^), (H) was presented as an example of 
a system for which the Smale horseshoe map can explic- 
itly be embedded as a subsystem into the flow along the 
homoclinic and heteroclinic orbits. This completed the 
proof of chaotic behavior of the studied system. 

Similar results hold for geodesic motion in arbitrary 
non-homogeneous pp-waves with the structural function 
of the form / = C(", where n > 4, i.e. for a general 
n-saddle potential (^ . It was shown in j2^ , that the 
decomposition into isolating cells Rj, j — 1,2, ■■■ , n, each 
containing exactly one of the basic unstable periodic so- 
lutions Ilj, is analogous to the case n = 3. Subsequently, 
the orbits Ilj were proven to be hyperbolic |2^ and the 
existence of homoclinic and heteroclinic orbits was es- 



tablished iQ. Again, given any biinfinite sequence, un- 
countably many orbits can be found which pass from one 
block containing Ilj to the other in the specified order. 

In order to independently support these arguments for 
chaotic behavior of geodesies in non-homogeneous pp- 
waves we investigate the motion also by a fractal method. 
Complementary to the analysis described above, we con- 
centrate on unbounded geodesies and we do not restrict 
to the same energy manifold E = const. 

Chaos is usually indicated by a sensitive dependence of 
the evolution on the choice of initial conditions. The co- 
ordinate independent fractal method (see for example , 
05 Bi 10) starts with a definition of different asymp- 
totic outcomes (given here by 'types of ends' of all pos- 
sible trajectories). A set of initial conditions is evolved 
numerically until one of the outcome states is reached. 
Chaos is uncovered if the basin boundaries which sepa- 
rate initial conditions leading to different outcomes are 
fractal. Such fractal partitions are the result of chaotic 
dynamics. As we shall now demonstrate, we observe ex- 
actly these structures in the studied system. 

We integrate numerically the equations of motion given 
by (H), (^. The initial conditions are chosen such that 
the geodesies start at r = from a unit circle in the 
{x, j/)-plane (due to the homogeneity of the n-saddle po- 
tential all other geodesies can simply be obtained by a 
suitable rescaling) . It is natural to parametrize the initial 
positions by an angle S [— tt, tt) such that x(0) = cos (j), 
2/(0) = sincf). In Fig. 1 we present typical trajectories 
of geodesies for n = 3,4,5, when i(0) = = 2/(0). We 
observe that each unbounded geodesic escapes to infin- 
ity (where for n > 3 the curvature singularity is located) 
only along one of the n channels in the potential. The 
axes of these outcome channels are given in polar coordi- 
nates by the condition cos(n0j) = —1, j = I, ■ ■ ■ ,n, and 
represent radial lines 'of steepest descent' since V — oo 
as p ^ oo most rapidly along them. (For non-zero ini- 
tial velocities more geodesies prefer one of the channels 
but the character of motion does not change significantly 

In fact, any unbounded geodesic oscillates around the 
radial axis (f>j — {2j — l)7r/n of the corresponding j-th 
outcome channel. Introducing A(/>j(r) — 4>{t) — (j)j we 
find asymptotically that p « [( § - 1 ) \/^!72 (r^ -r)] 2/(2-") 
as /9 ^ 00 and 

A^j{t) « (r, - t)"(Acos [61n(r, - r)] 

+Bsin[61n(T, -r)]) , (9) 

where a = i(n-|-2)/(n-2), b = ^y/7n'^ - An- 4/(n-2) 
and A, B are constants. As the geodesies approach the 
singularity at p = 00, r — > Ts, frequency of their oscil- 
lations around (j)j grows to infinity while the amplitude 
of oscillations tends to zero. We may call this effect a 
'focusation'. 
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FIG. 1. Geodesies starting from a unit circle escape to in- 
finity only along one of the n channels. 

Let us return back to our observation that all un- 
bounded geodesies approach infinity through only n dis- 
tinct channels. These represent possible outcomes of our 
system and we assign them symbol j which takes one of 
the corresponding values, j G {1,2,..., n} (thus, for ex- 
ample, j — 1 means that the geodesic approach infinity at 
p = oo through the first channel with the axis t^i = ir/n 
as T — > Ts > 0). From Fig. 1 we observe that in certain 
regions the function depends very sensitively on ini- 
tial position given by (j). We calculated numerically 
for n = 3,4, 5 — the results are shown in Fig. 2. Also, 
in the same diagrams we plot the function Ts(0) which 
takes the value of the parameter r when the singularity 
at jO = (X) is reached by a given geodesic. 



FIG. 2. The functions j{(f)) and Ts{(f>) indicate that basin 
boundaries separating different outcomes are fractal. 

Clearly, the boundaries between the outcomes appear 
to be fractal which can be confirmed on the enlarged de- 
tail of the image and the enlarged detail of the detail etc. 
In Fig. 3 we show such zooming in of the fractal inter- 
val localized around the value w for n = 3 (there 
are two symmetric fractal intervals in this case around 
(f) K, ±|7r) up to the sixth level. At each level the struc- 
ture has the same property, namely that between two 
larger connected sets of geodesies with outcome channels 
ji and j2 7^ ji there is always a smaller connected set of 
geodesies with outcome channel ja such that js ^ ji and 
js 7^ 32 ■ Similarly as in [Q, |]lO|, the structure of initial 
conditions resembles three mixed Cantor sets, and this 
fact is a manifestation of chaos. 
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FIG. 3. The fractal structure described by j{(f)) and Ts{<j)) 
is clearly confirmed here by zooming in the interval around 
(j) ~ for n = 3. 

The above structure of j (0) has its counterpart in the 
fractal structure of Ts{4>), see Fig. 2 and Fig. 3. We ob- 
serve that the value of goes to infinity on each dis- 
continuity of J (</»), i.e., on any fractal basin boundary 
between the different outcomes. There is an infinite set 
of peaks corresponding to chaotic bounded orbits which 
never 'decide' on a particular outcome, and so never es- 
cape to infinity. The value of Ts also increases in non- 
chaotic regions of as one zooms in the higher levels of 
the fractal. This is natural since these levels are given by 
geodesies which undergo 'more bounces' on the potential 
walls before falling into one of the outcome channels so 
that their values of Tg are greater. 

We demonstrated chaotic behavior of geodesies in non- 
homogeneous pp-waves by invariant analytic and numer- 
ical methods. As far as we know, this is the first ex- 
plicit demonstration of chaos in exact radiative space- 
times (chaotic interaction of particles with linearized 
gravitational waves on given backgrounds has already 
been studied in [|§ - Since pp-wave solu- 

tions are the simplest gravitational waves it would be an 
interesting task to search for a chaotic motion in other 
radiative space-times. 
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